Info<< "Reading field p\n" << endl;


volScalarField p
(
    IOobject
    (
        "p",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

Info<< "Reading field U\n" << endl;
volVectorField U
(
    IOobject
    (
        "U",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

Info<<"Reading fields T_1\n"<<endl;
volScalarField T_1
(
    IOobject
   (
    "T_1",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );

Info<<"Reading fields T_0_71\n"<<endl;
volScalarField T_0_71
(
    IOobject
   (
    "T_0_71",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );

Info<<"Reading fields T_0_25\n"<<endl;
volScalarField T_0_25
(
    IOobject
   (
    "T_0_25",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );
Info<<"Reading fields T_0_125\n"<<endl;
volScalarField T_0_125
(
    IOobject
   (
    "T_0_125",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );
Info<<"Reading fields T_0_05\n"<<endl;
volScalarField T_0_05
(
    IOobject
   (
    "T_0_05",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );
Info<<"Reading fields T_0_025\n"<<endl;
volScalarField T_0_025
(
    IOobject
   (
    "T_0_025",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );
Info<<"Reading fields T_0_005\n"<<endl;
volScalarField T_0_005
(
    IOobject
   (
    "T_0_005",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );
Info<<"Reading fields T_0_001\n"<<endl;
volScalarField T_0_001
(
    IOobject
   (
    "T_0_001",
    runTime.timeName(),
    mesh,
    IOobject::MUST_READ,
    IOobject::AUTO_WRITE
   ),
    mesh
 );
//*****************************ADD
Info<< "Reading transportProperties\n" << endl;

IOdictionary transportProperties
(
    IOobject
    (
        "transportProperties",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ_IF_MODIFIED,
        IOobject::NO_WRITE
    )
);

/*
dimensionedScalar DT
(
 "DT",
 transportProperties.lookup("DT")
 );
//
*/


//**********************************
#include "createPhi.H"
dimensionedScalar nu
(
    "nu",
    dimViscosity,
    transportProperties.lookup("nu")
);
dimensionedScalar Pr_0_001 (transportProperties.lookup("Pr_0_001"));
dimensionedScalar Pr_0_005 (transportProperties.lookup("Pr_0_005"));
dimensionedScalar Pr_0_025 (transportProperties.lookup("Pr_0_025"));
dimensionedScalar Pr_0_05 (transportProperties.lookup("Pr_0_05"));
dimensionedScalar Pr_0_125 (transportProperties.lookup("Pr_0_125"));
dimensionedScalar Pr_0_25 (transportProperties.lookup("Pr_0_25"));
dimensionedScalar Pr_0_71 (transportProperties.lookup("Pr_0_71"));
dimensionedScalar Pr_1 (transportProperties.lookup("Pr_1"));
dimensionedScalar gamma_0_001 (transportProperties.lookup("gamma_0_001"));
dimensionedScalar gamma_0_005 (transportProperties.lookup("gamma_0_005"));
dimensionedScalar gamma_0_025 (transportProperties.lookup("gamma_0_025"));
dimensionedScalar gamma_0_05 (transportProperties.lookup("gamma_0_05"));
dimensionedScalar gamma_0_125 (transportProperties.lookup("gamma_0_125"));
dimensionedScalar gamma_0_25 (transportProperties.lookup("gamma_0_25"));
dimensionedScalar gamma_0_71 (transportProperties.lookup("gamma_0_71"));
dimensionedScalar gamma_1 (transportProperties.lookup("gamma_1"));


label pRefCell = 0;
scalar pRefValue = 0.0;
setRefCell(p, piso.dict(), pRefCell, pRefValue);
mesh.setFluxRequired(p.name());


singlePhaseTransportModel laminarTransport(U, phi);

autoPtr<incompressible::turbulenceModel> turbulence
(
    incompressible::turbulenceModel::New(U, phi, laminarTransport)
);

#include "createMRF.H"
#include "createFvOptions.H"
